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The  Angular  Impact  Distibution  of 
Charged  Particles  Attracted  to  a 
Charged  Cylindrical  Spacecraft 


1.  INTRODUCTION 

Spacecraft  in  geosynchronous  orbit  exhibit  charging  when  immersed  in  hot 
substorm  plasmas.  The  potential  that  a  surfnce  attains  depends  both  on  the 
material's  properties  and  on  the  plasma's  characteristics.  Secondary  emission 
of  electrons  by  electrons  and  ions,  as  well  as  backscattering.  determine  the 
floating  potential.  These  are  both  functions  of  the  angle  of  incidence  of  particles 
on  the  surface. 

In  this  work,  the  angular  distribution  of  incident  particles  is  obtained  for 
attracted  and  repelled  particles  for  cylindrical  spacecraft.  Results  are  given 
for  the  plane  and  spherical  cases  as  well. 

For  isotropic  spacecraft  materials,  the  backscattering  coefficients,  the 
secondary  emission  coefficients,  and  the  shielding  efficiency  depend  on  impact 
energy  and  impact  zenith  angle,  but  are  independent  of  azimuth  angle.  According¬ 
ly,  only  the  distribution  of  impact  zenith  angles  is  important.  For  convenience, 
only  protons  will  be  considered,  the  extension  to  other  charged  particles  being 
obvious.  First  the  problem  will  be  precisely  defined,  then  a  simplified  model  of 
physical  reality  set  up,  and,  finally,  the  distribution  of  impact  angles  derived. 
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PHBCBDIW  JMOB  KAIK-NOT  flLMD 


2.  DKKI.MTION  OF  ^•KOI^I,KM 


Kind  tho  flux,  ,  of  pi  olons  slrikinn  an  i-lonn-nl  of  ai  i’a  i\ith  inipai  t 
kini'tii'  riKM'git'S  within  tho  band  K  ±  ^K'2  and  with  inipaot  /onilii  ancli  s  within 
tho  band  0  t 


2.  MODKI. 


riio  inodol  is  that  of  an  infinit  ivolv  long  I'vlindor  of  radius  It  .  Itovond  a 

I  > 

"shoath"  of  radius  l{^,  tho  plasma  is  isotropi<-  and  is  at  /oro  [lotontial.  I'ho 
potontial  within  tho  shoath  is  niodolod  by: 


ln(H/H^) 

o  ln(R  /R  ) 
o  s 


(1) 


I'ho  subsoripts  "o"  and  "s"  will  bo  used  to  mdioato  valuos  of  variables  at  R  and 
R  ..  No  subsoript  will  bo  used  for  gonin-al  valuos  of  H. 

In  tho  model,  tho  differential  flux  f(K,  <l>)  piotons  cm  “  see  stor  oV  is 
a  function  of  kini'tio  and  pot<'ntial  onorgios  only  —  oxi'opl  in  tho  "shadow"  of  the 
satellite  whore  it  is  zero.  Both  potontial  (4)  and  kinotn-  (K)  onoi  gy  are  expressed 
in  elootron  volts.  Any  impact  angle  and  energy  oorrosponding  to  a  trajoi'lory 
originating  elsewhere  on  tho  satellite  is  oonsidorod  to  bo  "shadowed." 


4.  DKRIV  ATION  OK  IMPACT  /KNITIl  ANGLK  OISTRIBl  TION 

The  impart  proton  flux  is  th<'  product  of  four  factors,  namely:  energy  band 
width,  unshadowed  solid  angle,  projected  area,  and  differential  flux.  This  may 
be  w  ritten  as 

i  -  (AK)(A0  <(i  sin  0  )(^A  eos  0  )(f(I-;  *  ))  (2) 

•a  z  o  o  o’  o 

where  =  total  unshadowed  azimuth  angle.  The  solid  angle  is  that  which  would 

be  subtended  at  the  earth's  center  by  the  area  enclosed  by  latitudes  0  i  AS  2  and 

■  ■  o 

longitudes  0  i  azimuth  angle  will  be  taken  as  zero  perpendicular  to  the 

cylinder's  axis.  If  shadowing  occurs,  it  will  be  at  the  larger  angles. 

Flquation  (2)  may  be  rewritten  as 

j  (AE  A0AA  f(E  ,  *  )l[sin  0  cos  0  )  (4(*  |  (3) 

■'o  o  o  o  o  '  *^m 


fi 


wheie  -  maximum  unshadowed  azimuth  angle,  di  being  measut-ed  So  as  alwavs 
to  lie  in  the  fii  st  quadrant.  The  faetor  4  accounts  for  the  four  quadi  anis.  A 
further  rewriting  yields 

i  =  K  di  sin  2  0 
■a  in  o 

K  =  4^  Ad  AA  f(K  ,  $  )  (4) 

0  -c  i  <  71 '2  . 

ni 

Note:  sin  2  0  s  2  sin  0  cos  0 
o  o  o 

Equation  (4)  is  very  general.  The  geometry  of  the  spacecraft's  surface  and  field 
enter  into  the  equation  only  through  d',.|.|,  the  maximum  unshadowed  azimuth  angle. 

5.  SPECIAL  CASES 

For  positive  spacecraft  potentials,  the  reverse  tra.iectories  curve  away  from 
the  spacecraft;  therefore,  there  is  no  shadowing  and 

<4  ■-  )r/2  if  $  ?  0  .  (a) 

From  Eq.  (4),  the  angular  distribution  is: 

i  =  (ir/2)  K  sin  2  0  . 

■  o  ' 

The  other  special  case  is  that  of  looking  for  protons  of  impossibly  low  energy. 
Here  we  may  write 

<4  -  0  if  $  <  -E  ,  (d) 

m  o  o 

which  has  meaning  only  for  negatively  charged  spacecraft. 


6.  MAXIMUM  U^SHADOWED  AZIMUTH  ANCLE 

The  simplest  cases  have  been  dealt  with  in  the  preceding  section.  This 
section  will  consider  cases  defined  by: 


(7) 


-E  <  $  <0 

l)  o 


The  model  has  zero  axial  field.  Therefore,  the  kinetic  energy  associated  with 
the  axial  velocity  component  is  independent  of  radial  position  and  may  be  written: 

E  -  E  sin^  (j)  sin“  6  .  (8) 

a  o  o  o 

where  E^  is  the  impact  energy  and  and  6^  are  the  impact  angles.  This  rela¬ 
tionship  follows  from  the  breakdown  of  the  impact  velocity  into  four  components, 
as  shown  in  Figure  1.  The  model  has  a  radial  field  about  the  cylinder's  axis. 
Thus  the  tangential  (i  to  radial  and  axial  directions)  kinetic  energy  is  subject  to 
conservation  of  angular  momentum  and  is  given  by: 

E,(R)  =  E  (R  /R)^  cos^  <t>  sin^  0  .  (9) 

t  o  o  o  o 

The  radial  kinetic  energy  is  simply: 


E^  =  E  cos^  0  .  (10) 

The  total  proton  energy  is  independent  of  the  radius,  and  beyot:d  the  sheath  it 
equals  kinetic  energy.  The  total  energy  is  thus: 

E  =  E+  E+E,  +  $  ,  $<0  .  (11) 

I'  a  t 

At  the  sheath's  edge  Eq.  (11)  becomes: 

E  =  E  cos^  0  +  E  sin^  d  sin^  0  +  (R  /R  )^  K  cos^  d  sin^  0  (12) 

s  so  o  o  o'  s  o  o  o 

and  at  the  satellite's  surface  it  becomes 

E  =  E  cos^  0  +  E  sin^  d  sin^  0  +  E  cos^  d  sin^  0  +  $  (13a! 

o  oc.  o  oo  o  oo 


E,^  =  E  +  ® 
s  o 


4  <  0  . 


The  limiting  case  corresponds  to  a  vanishing  radial  energy  at  the  sheath's  edge  or 


cos  0=0 
s 


I 
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Equations  (4>  and  (lii)  together  give  the  zenith  angle  (0)  distribution  of  protons 

striking  the  surface  c-f  a  cylindrical  satellite. 

The  angular  distribution  derivation  given  has  considered  only  the  potential  at 

the  inner  (R  =  R  )  and  outer  (R  =  R^)  edges,  not  the  potential  within  the  sheath 

(R  <  R  ).  It  may  be  shown  that  the  potential  within  the  sheath  is  noi  controlling, 

'  1 
provided  that  at  all  points  within  the  sheath  it  is  equal  to  oi-  more  negative  than 

that  of  the  inverse  square  potential  function: 


(R  /R)^- 1 
s _ 

(R  ,/R  )“- 1 


<1>  <  0 
o 


(17) 


This  function  is  adjusted  to  match  the  model  at  the  sheath's  edge.  The  logarith¬ 
mic  function  of  the  model  satisfies  this  condition.  The  proof  is  similar  to  that 
given  for  the  spherical  case  for  objects  in  a  gravitational  well.  The  proof  is  only 
valid  if  the  potential  is  a  function  of  radial  position  alone. 

The  distribution  of  impact  zenith  angles  according  to  Eqs.  (4)  and  (lii)  is 
shown  in  Figure  2  for  several  cases.  The  deviations  from  the  symmetrical 
sin  (2  0)  distribution  are  due  entirely  to  shadowing.  The  shadowing  ;s  equivalent 
to  a  "horizon"  at  less  than  90”  below  the  zenith  due  to  refraction  of  tlie  protons. 


0“  10"  20*  30*  40*  50*  60°  70°  80°  90° 


ZENITH  ANGLE 

Figure  2.  Distribution  of  Zenith  Impact  Angles  for  Protons  on  a  Cylindrical 
Satellite:  (a)  Thin-Sheath  C  ase 
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T.  DISCI  SSI()> 


rt'.i'  •hi'iii  y  civi'ii  hiTf'  is  I'XMcl  iii  !li;il  il  is  i  i 'inpli  ' .-K  i  i .:i,si.s' •  ifh  Mu- 
i;:,ul,  l.  I'hi'  ir.odi‘1,  til i\v I’Vc'i',  1  ^•p!■^■sotlt.s  a  pliysii  al  iiapi issi la  1 1 1 y  in  ilial  lla  '■li'<  - 
M  ia  fii'lrl  IS  disi'oiilinunus  at  flic  sui  facc  cf  llic  slicatli  and  an  appn ixnnat  i ■  m  in 
tlia'  end  effects  at  e  taken  to  be  vanisliinely  small,  (inly  if  the  tu’afest  end  is  at  a 
distance  ceev  laege  cempafed  to  the  radius  uill  its  effect  become  eanisliiiigly 
small.  file  theory  does,  however,  represent  :t  limit  ing  case  and,  as  such,  is 
useful. 


8.  SI  MM  VKV 


The  angular  distribution  of  impacting  repelled  panicles  is  sin  (2  '0,  as  given 
by  i:q.  (4),  since  there  is  no  shadowing  for  repelled  parti. des.  electrons  inci¬ 
dent  on  a  negatively  charged  surface  and  ions  on  a  positively  chargi'd  surface  nill 
im[)act  in  this  way.  In  order  to  ealculate  surfaei'  potentials,  taking  aci  ouni  'f 
sei  .iiidary  electron  emission,  one  must  know  the  distribution  of  impact  angles  a.s 
well  as  the  seco.idary  electron  yield  as  a  function  of  angle. 

The  angular  distribution  of  impacting  tittracied  t'.artieles  (electrons  on  a 
positive  sut'ifare,  and  piusitive  ions  on  a  negative  surface)  is  mote  complicated. 

The  nipacting  angular  distribution  for  a  cylindrieal  surface  for  attracted  parlicle.- 
has  been  derived  here,  as  shown  in  Figures  2  (a),  2  (b),  and  2  (i  ),  for  a  r.angi-  of 
energies  and  sheath  thickness. 

Till'  distribution  of  impact  angles  (attracted  particles)  for  a  plane  surface  as 

'111  same  .as  'he  ease  of  zero  sheath  thickness.  For  this  ease  the  distribution  is 

[iri iport ional  to  sin  (2  0)  for  0  less  than  0  ,  where 

m 

-  cos"^  (4)/E)^'“  , 

and  where  F  is  the  impact  energy.  If  0  is  greater  than  the  distribution  func¬ 

tion  IS  zero.  For  spherical  surfaces  there  is  no  shadowing,  so  the  distribution  :s 
again  proportional  to  sin  2  0  in  the  thick  sheath  limit. 
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